The Harmonic Measure of Diffusion-Limited Aggregates including Rare Events 
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We obtain the harmonic measure of diffusion-limited aggregate (DLA) clusters using a biased 
random- walk sampling technique which allows us to measure probabilities of random walkers hitting 
sections of clusters with unprecedented accuracy; our results include probabilities as small as 10 -80 . 
We find the multifractal D(q) spectrum including regions of small and negative q. Our algorithm 
allows us to obtain the harmonic measure for clusters more than an order of magnitude larger than 
those achieved using the method of iterative conformal maps, which is the previous best method. 
We find a phase transition in the singularity spectrum f(a) at a w 14 and also find a minimum q 
of D(q), q mm =0.9 ±0.05. 
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I. INTRODUCTION 

Diffusion-limited aggregation (DLA) is a stochastic 
model for irreversible growth which gives rise to frac- 
tal clusters [H, 0| , see Figures 1 1121 The growth process is 
defined by releasing a random walker far from the clus- 
ter and allowing it to diffuse until it sticks to the surface 
and becomes part of the cluster. Then another particle 
is released, and so forth. The probability of sticking at 
various points on the cluster, i.e. the distribution of the 
growth probability, is a function with very large varia- 
tions. It is the subject of this paper. 

Since the Laplace equation is equivalent to the steady- 
state diffusion equation, this probability distribution is 
proportional to the perpendicular electric field on the 
surface of a charged electrode with the shape of the clus- 
ter; in this context the probability is called the harmonic 
measure, and is defined for any surface. For fractal sur- 
faces, including that of DLA, the harmonic measure is 
usually multifractal [3[. For DLA the harmonic measure 
is of particular interest because of the connection with 
the growth probability. For other fractal surfaces this 
connection is lost. However, the measure is still of sub- 
stantial practical interest because its relationship with 
physical processes such as catalysis [4j. 

For many interesting equilibrium fractals the harmonic 
measure can be calculated using conformal field theory [5, 
@, 0] or Schramm-Loewner evolution (SLE) [§[ . There is 
no corresponding theory for DLA for which the measure 
must be found numerically. There are numerous studies 
in the literature of this quantity, for example 0, 
fl2j . This is a difficult problem because of the very large 
variation of the growth probability. As we will see the 
dynamic range of the function is of the order of 10 80 even 
for rather small clusters. This is far out of the range 
accessible to straightforward random walker sampling. 

In this paper, we use a biased random-walk sampling 
method. We can obtain extremely small growth proba- 
bilities and to accurately obtain the complete harmonic 
measure for DLA clusters of up to 10 6 particles. The 



method was previously used on percolation and Ising 
clusters [l]|. For those (equilibrium) systems, we found 
good agreement with analytic predictions for the har- 
monic measure [3, [l5| . 

The harmonic measure is usually characterized in 
terms of the generalized dimensions D(q). For integer 
q, D{q) corresponds to the fractal dimension of the q 
point correlation function. We define D(q) by partition- 
ing the external boundary of a DLA cluster into boxes of 
length I. The probability that a diffusing particle will hit 
the section of the perimeter contained in box i is denoted 
by pi. These probabilities define a "partition function" 
Zl(q) = J2i(Pi) 9 P3- ^ Ziil) can be written as a power 
law in the dimensionless ratio R/l, where R is the over- 
all size of the cluster, then the generalized dimension is 
given by: 



Zi(q) = (R/l)- T « = (R/l)- (q - 1)D{q) . 



(1) 



There are special values of D(q) including the box- 
counting dimension, D(0) f=a 1.71 [l7T |. For two dimen- 
sional clusters we always have D(l) = 1 [Lsl ]. 

Another quantity of interest is /(cv), which is called 
the singularity spectrum. This function is the Legendre 
transform of r(q): 



j,, \ dr dr 
= q— -t, a=—. 
aq aq 



(2) 



As is the case for D(q), some special values of /(a) are 
known: f(l) = 1 and the largest value of f(a) is equal 
to D(0). /(a) can have a phase transition, namely a 
maximum value of a, a max , for which f{a) is defined. 
There has been significant disagreement as to whether 
/(a) for DLA has a phase transition. This controversy 
is summarized in [ll[. 

The main difficulty in resolving this issue is that large 
a, or small q, corresponds to the smallest probabilities 
on the cluster. The straightforward method of obtain- 
ing the harmonic measure, sending large numbers of ran- 
dom walkers at the cluster, is only capable of measur- 
ing probabilities down to ~ 1CP 10 ; even clusters with 



only 1 , 000 particles have sections with growth proba- 
bilities significantly smaller than that. This issue was 
partially resolved by Jensen et al. [HI], wn0 used the 
method of iterated conformal maps (CM) HHHHll, to 
obtain significantly smaller probabilities. Their main re- 
sult was the determination of the harmonic measure of a 
single cluster of size 3 TO 4 , where they found probabilities 
down to 10~ 35 . This work [ll[ was a significant advance, 
though the CM results are not conclusive in giving the 
asymptotic results for DLA because the CM method is 
limited to small clusters, and it is known that some fea- 
tures of DLA have slow crossover to asymptotic scaling 
[22| . There are other, technical, questions about the CM 
method that we discuss below. 

Other groups have obtained the harmonic measure for 
on-lattice clusters using relaxation methods to solve the 
Laplace equation. Ball and Spivack [l(| grew DLA clus- 
ters, corrected for lattice anisotropy, up to 10 5 particles. 
They then solved Laplace's equation numerically to ob- 
tain the measure. Hanan et al. [12] measured the com- 
plete harmonic measure of DLA clusters using a related 
relaxation technique. In contrast to [ic| |. these authors 
first grew the cluster off-lattice, then forced it on lattice 
to solve for the measure. The simulations in [12] were 
also limited to small clusters of 6 ■ 10 4 particles. 



II. SIMULATION METHODS 

We grow our DLA clusters by the method that is now 
standard [2|, which includes speeding up the process by 
allowing the random walker particles to take large jumps. 
We store the cluster in a data structure which allows the 
calculation of the size of the jump to be performed in 
0(log(n)) time, where n is the cluster size. These meth- 
ods allow us to grow clusters in 0(nlog(n)) time, a big 
improvement over the CM method, which is 0(n 2 ). In 
the CM method the harmonic measure is available at each 
step. In our case, we use a biased random walker method 
(the signpost method [HI]) to obtain the harmonic mea- 
sure once the cluster has finished growing. 

The signpost method consists of two iterated steps: 
a sampling step and a measurement step. On the first 
step, a large number of random walker probe particles 
(N), each with weight 1/N, are released far from the 
cluster and diffuse until they hit the cluster. This allows 
us to determine which areas of the cluster are poorly 
sampled. Next, we place line segments (signposts) block- 
ing off all regions the cluster that have sites that are hit 
by fewer than some percentage of the probe particles, 
say 10%. In the measurement step, we release N more 
probe particles far from the cluster and allow them to hit 
the cluster and signposts. The probe particles that hit 
the cluster have their weight permanently added to the 
perimeter site probability distribution. The locations on 
the signposts where the probe particles hit in the first 
step are used as the initial location of the N probe parti- 
cles released in the probe step of the second iteration. To 
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FIG. 1: A branch of a DLA cluster with an external border 
representing the Minkowski cover of the cluster. The particles 
filled grey (blue online) in the center are accessible to random 
walkers. The solid black particles can never be reached. Inside 
the gray box is a narrow neck which marks a low probability 
pathway for random walkers. 



conserve probability, each probe particle released in the 
second iteration has weight p/N, where p is the fraction 
of probe particles that absorbed on the signpost in the 
first iteration. The probe particles released in the probe 
step of the second iteration help determine which sections 
of the cluster are still poorly sampled. More signposts 
are added to block off the still poorly sampled regions 
and then the probe particles for the measurement step 
are released. This process is repeated until the growth 
probability of all sites has been measured. For a more 
detailed description of the algorithm, see [T^] ■ 

This method is similar to a rare event method in chem- 
ical physics that uses 'milestones' [23j]. The main differ- 
ence between the our method and that of [23| is we choose 
the locations of the signposts/milestones dynamically 
and that we do not need to know an a priori distribution 
for the random walkers along the milestones. See also 

This signpost method allowed us to measure probabil- 
ities down to 10 -300 for percolation and Ising clusters. 
For DLA we have measured probabilities down to 10 -80 . 
Performing the signpost algorithm on DLA clusters is 
more complex than the percolation and Ising cases. DLA 
clusters are grown off-lattice, which means that some sec- 
tions of the exterior of the cluster are almost completely 
blocked by two branches of the cluster nearly meeting, 
making a narrow passage slightly larger than the diame- 
ter of a probe particle; see Figure [1] The probability of 
a probe particle diffusing through some of these passages 
without touching the cluster is smaller than 10~ 8 . These 
passages arc treated differently, but in a way consistent 
with the signpost algorithm. Specifically, we slowly move 
signposts perpendicular to the passage inward over sev- 
eral iterations until the probe particles can reach other 
sections of perimeter. In other words, we allow narrow 
passages to have closely spaced signpost lines so that 
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FIG. 2: The Harmonic Measure for a cluster with 10 4 parti- 
cles, the lighter the color the smaller the measure. The size 
of the particles is doubled to represent the cover of the DLA 
cluster. The smallest measure in the cluster is ft! 10 -49 . Sites 
that cannot be reached are marked grey. 



proper sampling can be achieved. 

One minor difference between the signposting we use 
in this paper and the one used previously is that we now 
reduce the threshold for blocking off sections of the clus- 
ter as a function of the total number of probe particles 
that hit signposts in the previous iteration. When more 
probe particles hit the signposts, we move the signposts 
much deeper the next time. Previously, we reduced the 
threshold a fixed amount each iteration. We found that 
this dynamic threshold adjustment gave us more consis- 
tent particle saturation on the signposts, which in turn 
decreased the rate at which the error grew from step to 
step. 

Before we can apply the signpost method to a DLA 
cluster, we first must find the perimeter of the cluster. 
More precisely, we must find all sites that are accessible 
to the probe particles. First, we take a ball the size of 
a probe particle and roll it clockwise around the cluster 
particle furthest from the center of the cluster until the 
ball touches a second cluster particle. After that, the ball 
is rolled clockwise about the second cluster particle until 
it touches a third particle. This process is repeated until 
the ball returns to its initial location. Note that a sin- 
gle cluster particle can be visited more than once by the 
ball. This process finds something akin to the Minkowski 
cover of the cluster, Figure [TJ We found that on average 
only 80% of the particles in cluster are accessible to ran- 
dom walkers. This means that 20% of the particles have 
a measure of exactly zero and these regions of the cluster 
will never grow. We found the 80% accessibility to be 
constant over a range of large cluster sizes, which shows 
that the accessible perimeter has the same fractal dimen- 
sion of the complete perimeter and the cluster itself, in 
contrast to percolation where the corresponding accessi- 
ble perimeter has a smaller fractal dimension than the 
complete perimeter [l4], HI] ■ 

Our perimeter accessibility results agree qualitatively 
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FIG. 3: f(a) vs. a for seven different system sizes with error 
bars. Note that the spectra appear to be converging to some 
asymptotic spectrum. 
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FIG. 4: /(a) vs. a for n = 10 6 and the asymptotic estimate, 
dotted, with error bars estimated from the data in Figure [3] 
Note that the asymptotic estimate terminates at a = 11 only 
because there were too few system sizes to extrapolate for 
larger a. We believe the phase-transition in f(a) occurs at 
a w 14. 



with other work [26| , which looked at cluster particle ac- 
cessibility as a function of probe particle size. However, 
the agreement is not precise because the authors of [26[ 
measured the accessible cluster using 10 5 random walk- 
ers, which are extremely unlikely to hit the low measure 
sections of the perimeter. 

We grew DLA clusters of various sizes: 10 3 , 3-10 3 , 10 4 , 
3 • 10 4 , 10 5 , 3 • 10 5 , and 10 6 particles. For each DLA clus- 
ter grown, we obtained the harmonic measure using the 
signpost algorithm. Figure [5] shows the harmonic mea- 
sure for a cluster with 10 4 particles. The different cluster 
sizes required a different number of random walkers per 
iteration, 10 6 , 10 6 , 5 • 10 6 , 10 7 , 2.5 • 10 7 , 10 8 , and 2.5 • 10* 
for 10 3 , 3-10 3 , 10 4 , 3-10 4 , 10 5 , 3-10 5 , and 10 6 sized clus- 
ters, respectively. The number of random walkers needed 
was estimated by determining the number of walkers re- 
quired to get at least 10 4 random walkers absorbed on 
each signpost for every iteration. We believe this is a 
conservative criterion. 
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FIG. 5: The slope of the power-law fit to the probability 
distribution at various points. The slope at each probability 
is averaged over about an order of magnitude in probability. 



III. RESULTS 

We use the method described above to obtain D(q). 
First, we take the space that contains a cluster and sec- 
tion it into boxes of size / and then measure Zi(q). We 
do this measurement of Z\ (q) for various values of I for 
a given q. Next, we calculate the slope of the function 
\nZi(q) versus InZ; this is r(q), which when divided by 
(q— 1), gives D(q). The fit is performed over the range of 
I for which the log-log plot is linear. This range is about 
one order of magnitude for the smallest system size and 
larger than one order of magnitude for larger systems. 
With this set of D(qYs for individual clusters of various 
sizes, we can calculate the average values of D(q) for var- 
ious sizes. We found that our results for large systems 
are close to the known values for -D(O) and D(l), 1.66 
and 0.99 respectively. 

Using the results of D{q) for individual clusters, we can 
Legendre transform the results to obtain f(a) for each 
cluster; see figure [3] For a range of a's, we estimated 
the asymptotic value of f(a) using finite-size scaling tech- 
niques in n, so that the correction to f{a) is of the form 
where j3 is a crossover exponent. The asymptotic 
values were determined by minimizing the residual of the 
power-law fit, see Fig. @] We found best fit for j3 was 0.4 
for the entire range of a. Note that exponents of 1/3 
and 1/2 are also consistent with the data. This means 
we are consistent with [22J where a crossover exponent of 
1/3 was found (for different quantities). The asymptotic 
f(a) values are consistent with special known values /(l) 
and the maximum /, measured to be 1.00 and 1.71 re- 
spectively. We believe the asymptotic f(a) calculated is 
the true f(a) for DLA. The last a for which f(a) is de- 
fined is more difficult to estimate. From visual inspection 
the asymptotic point of the phase transition appears to 
be about a rj 14. This is significantly smaller then the 
value found by the authors of [ll|, a ~ 18. 

We were also able to obtain a histogram of the growth 
probabilities for every system size. The bins were sized 



logarithmically, to allow for a power-law fit to the results. 
Figure O shows the slope of the power-law fit to the prob- 
ability distribution. The slope is fairly consistently —0.9 
which corresponds to a smallest q for which D(q) is de- 
fined being —0.1. These values agree moderately well 

with [n|. 



IV. CONCLUSIONS 

In this paper we applied a rare-event technique to 
obtain the complete harmonic measure of DLA clus- 
ters. We found that probability distribution is consis- 
tent with a power-law exponent of -0.9. We also found 
a slow crossover to infinite-size cluster behavior in /(a), 
in agreement with previous work (l2l . [22| . We believe 
that our extrapolated f{a) is a very good approximation 
to f(a) for infinite-sized DLA. We found a phase transi- 
tion in f(a) at a ks 14. This maximum a is related to 
the opening angle of the branches near the seed point of 
the cluster. The area around the seed point should have 
the lowest measure, so the angles in that region, 4 > mn x 
are related to the largest alpha by a max ~ l/4> m ax, \XM- 
By determining the exact relationship between a max and 
4>max, Hanan and Heffernan |12j determined the asymp- 
totic a max as a max « 15 using results from Mandelbrot 
et al. [27| for the asymptotic estimate of <f> ma x- This is 
in satisfactory agreement with our results for a max . 

Our results differ significantly from those obtained by 
the CM method [ll[ in several ways. First, we find a sig- 
nificantly smaller value for a max . Second, we find that 
finite-size effects are still noticeable on clusters with 10 6 
particles. The authors of [ll[ found no finite-size effects 
at their largest system size, 3 TO 4 particles. This is incon- 
sistent with our findings. Lastly, we find that the smallest 
probabilities found on clusters of size 3 • 10 4 are signifi- 
cantly smaller, about 5 to 10 orders of magnitude, than 
reported for CM clusters. We do find good agreement 
for small and moderate values of a, which corresponds 
to region of the spectrum which is easily measured by 
random walker sampling. This explains why the differ- 
ence between CM and standard DLA clusters was not 
seen earlier. 

Assuming that the signpost method and the CM 
method are both successful at obtaining the measure for 
their respective clusters, then the only explanation for 
the discrepancy is that CM clusters are not the same as 
DLA clusters grown using particles. Superficially, CM 
clusters appear to be the cover of DLA clusters. If this 
were the case, then both methods would obtain the same 
measure because the measure for a probe particle the 
same size as a cluster particle hitting a standard cluster 
is exactly the same as the measure for a point-sized par- 
ticle hitting the cover of the same cluster. The heart of 
the issue may be the shape and size of the 'bumps' added 
to the CM clusters during each step. These bumps are 
designed to have a semicircular shape, and to be of fixed 
size, but they can distort as noted in (28l . [29| . It is im- 
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portant to check that the bumps are, in fact of fixed size, Science Foundation through DMS-0553487. 

and resize them if necessary. It is not clear that this 

was done in Even if this correction were made, the 

shape of the bumps can be very distorted deep inside the 

cluster. 
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